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ABSTRACT 

The Energy Density Functional theory is one of the most used methods developed in nuclear structure. It is based 
on the assumption that the energy of the ground state is a functional only of the density profile. The method is 
extremely successful within the effective force approach, noticeably the Skyrme or Gogny forces, in reproducing the 
nuclear binding energies and other bulk properties along the whole mass table. Although the Density Functional is in 
this case represented formally as the Hartree-Fock mean field of an effective force, the corresponding single particle 
states in general do not reproduce the phenomenology particularly well. To overcome this difficulty, a strategy has 
been developed where the effective force is adjusted to reproduce directly the single particle energies, trying to keep the 
ground state energy sufficiently well reproduced. An alternative route, that has been developed along several years, for 
solving this problem is to introduce the mean field fluctuations, as represented by the collective vibrations of the nuclear 
system, and their influence on the single particle dynamics and structure. This is the basis of the particle-vibration 
coupling model. In this paper we present a formal theory of the particle-vibration coupling model based on the Green’ 
s function method. The theory extends to realistic effective forces the macroscopic particle-vibration coupling models 
and the (microscopic) Nuclear Field Theory. It is formalized within the functional derivative approach to many-body 
theory. An expansion in diagrams is devised for the single particle self-energy and the phonon propagator. Critical 
aspects of the particle-vibration coupling model are analysed in general. Applications at the lowest order of the 
expansion are presented and discussed. 

PACS : 2I.60.Jz, 0.5.30.Fk 


I. INTRODUCTION 

One of the most striking features of nuclei is the validity of the shell model where the nucleons move independently 
in a smooth single-particle (s.p.) potential.This potential is static, but non-local in space, which implies that the 
mean field is velocity (momentum) dependent. The fluctuations of this average potential leads to collective modes,in 
particular surface vibrations (phonons) [T|. Taking into account the coupling of these modes to the single particle 
motion, as in the particle-vibration coupling (PVC) model, the shell model acquires a dynamical content. Thus, the 
average potential becomes also non local in time, being characterized by an energy (frequencyIndependence. The 
dynamics of the shell model (the expression ’’dynamics of the shell model” is after G.E. Brown [^) affects different 
nuclear properties as fragmentation and related spectroscopic factors of s.p. states, their density (proportional to the 
effective mass m* near the Fermi energy), the s.p. spreading widths, the imaginary component of the optical potential. 
A unified description of the s.p. motion at positive and negative energy eventually emerges Q. The importance of 
the dynamical coupling was first demonstrated in the case of electrons in metals and of normal liquid ^He. These 
results inspired the work of Bertsch and Kuo @ on the enhancement of the effective mass near the Fermi energy in 
finite nuclei. The review article in ref. Q contains the many results obtained in the ’80. A recent introduction to the 
subject can be found in ref. Extensions of the PVG approach to the continuum @ and to the scattering processes 
0 have been presented in the literature. 
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Most of the calculations done in the past are based on purely phenomenological inputs. In Tables 4.3a and 4.3b 
of Ref. 0, an extensive review of the results obtained in by nine groups in the period 1968-1983 can be 

found. Rather different frameworks had been adopted, no one being fully self-consistent: s.p. potentials range from 
harmonic oscillator (HO) to Woods-Saxon (WS) or Hartree-Fock (HF) with Skyrme forces; residual interactions 
at the particle-vibration vertex are either multipole-multipole forces, or forces of Landau-Migdal type, or Skyrme 
forces but with velocity-dependent terms dropped, or even G-matrix interactions. Consequently, although there is 
qualitative agreement among several calculations, it is rather hard to assess seriously the quantitative impact of all the 
approximations done and compare in detail with EDF-based calculations with so called spectroscopic accuracy and 
with the recent relativistic calculations These last one are fully self-consistent , because the whole interaction 

is considered both to build the phonons within the relativistic random phase approximation (RRPA) framework and 
to construct the PVC vertex. Within the framework of non-relativistic effective interactions, a step in this direction 
has been done very recently in ref. [ll|) [13 with calculations of ^°Ca and ^°®Pb s.p. self-energies based on a fully 
self-consistent random phase approximation (RPA) of the vibrations, and without dropping any term in the coupling 
vertex. 

In this paper, along the scheme presented in ref. [l3, we derive coupled integral equations for the single particle 
self-energy and Green’ s function, the phonon propagator and the vertex function. These equations form a closed 
system and we propose a solution by iteration. The terms so obtained by iteration can be represented by diagrams, 
that automatically contain the correct symmetry factors, as well as the Pauli corrections demanded by antisymmetry 
between the explicit particle lines and the ones implicitly contained in the phonon propagator. We propose a solution 
for the single-particle self-energy at the lowest order of approximation, with numerical applications for both density- 
independent and density-dependent forces. 

The outline of the paper is as follows. In Section [H] we present the basic formalism, starting from the construction 
of the hamiltonian in second quantization, from which we derive the equation of motion for the single particle Green’s 
function and the integral equations for the phonon propagator and the vertex function. For simplicity in this section 
the treatment is restricted to density independent forces. In Section iHll the expansion in diagrams is devised following 
an iterative procedure. Here some peculiarities of the formalism, which are consequences of the microscopic structure 
of the phonons, are discussed in some detail. Section m is devoted to the case of a density dependent force. We 
show that, at least for the applications at lowest order that are discussed in this work, this case can be reduced to 
the density independent case if the phonons are assumed to be of small amplitude. In Section |V] we present the 
application of the formalism to the nucleus ^'^Ca, and compare the results for a density independent and a density 
dependent Skyrme force. Conclusions are drawn in Section IVII 


II. THE FORMALISM 


In the general Energy Density Functional (EDF) approach it is assumed that the energy of the ground state can 
be obtained by the minimization of a functional that depends only on the single particle density profile. The method 
is based on the Hohenberg-Kohn theorem that proofs the existence of such a functional, although it does not 
provide any procedure to construct it. The theorem has been recently extended [IM3 to self-bound systems like 
nuclei. In nuclear physics the most widely used scheme to formulate the functional is based on the introduction of 
effective forces, noticeably Skyrme or Gogny forces [13. These forces are phenomenological and are assumed to be used 
in the mean field scheme and to incorporate part of the effect of correlations through phenomenological parameters 
to be fitted to reproduce the binding energy of nuclei, as well as the trend of radii and deformations, throughout 
the mass table. Other functionals, e.g. the recent BCPM [l9l - l?ij| . contains more microscopic input, but they still 
include phenomenological terms to be fitted. Not only masses, but also other simple functionals of the density like 
radii etc. can be well reproduced. These functionals have also shortcomings, and in this work we focus on the fact 
that observables like the single-particle states and their fragmentation cannot be accounted for. The approach has 
been extremely successful and an enormous literature exists on the subject. Although in what follows we do not need 
the particular form of the functional, let us consider the case of a functional that can be formally written as the 
Hartree-Fock mean value of an effective force /, 

E = T + -'^Pikfu-,kj{p)Pij = T + V{p), (1) 

iklj 


where T is the kinetic energy, pik are the matrix elements of the density, and ijk . are labels for generic single 

particle states. The effective force / is in general a functional of the density matrix pik- The matrix elements of the 
effective force are antisymmetrized 


rA 

J il\kj 


< il\f\kj > - < il\f\jk > . 


( 2 ) 
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The density matrix can be written as the expectation value 

Pik =<'ipt Vfe > (3) 

where are the single particle creation and annihilation operators. The density dependence can originate from 

three-body (or higher) forces, if one replaces a pair by its ground state average, Eq. m, thus reducing the 

three-body force to density dependent two-body force. However in general the force contains non-integer power of 
the density, and it cannot be derived from a many-body force. Notice that the antisymmetrization in Eq. (I2])is not 
performed with respect to the orbitals of the density matrix that appears implicitly in the effective force. This is 
a general problem for the effective density dependent force, discussed by several authors [12, [12|, but we follow the 
common procedure to use expression O- We will develop the formalism assuming a density independent force like 
SV [15 . Later we will discuss how to implement the approach to the case of a density dependent force and we will 
discuss the problems that arise in that case. There are currently several attempts to fit functionals that are density 
independent, with the goal of avoiding problems that have been identified in the multi-reference DFT [25j |. 

The dynamics in the single particle degree of freedom can be described by the one-body Green’s function and the 
corresponding equations of motion. To this purpose one needs a Hamiltonian formulation. We will assume that the 
independent single particle motion is determined by the single particle potential, while the coupling with the density 
fluctuations is determined by the same interaction appearing in the functional of Eq. o, with / independent of 
density. In other words the Hamiltonian will contain a one-body term, which includes the single particle potential, 
and a residual two-body interaction which describes the fluctuation dynamics around the ground state. Therefore we 
will assume the following form of the Hamiltonian 

H = < k l \ f \m n >A N (ipl'ipj'ijjniprn) (4) 

klmn 

The second quantization form of the Hamiltonian is used in order to incorporate in a systematic way the Pauli 
principle in the equation of motion. The symbol indicates the normal product of the operators with respect 

to the independent particle ground state. This form ensures that the two-body interaction has vanishing mean value 
in this ground state. The term Hq is the one-body part of the Hamiltonian, that in general can be written 

Hq = ( 5 ) 

ik 

Here Cik is the single particle energy matrix, obtained from the Hartree-Fock approximation. This one-body part 
describes the single particle motion in the independent particle limit, i.e. when the two-body interaction term is 
neglected. The form of Eq. ([4]) for the hamiltonian in the equation of motions poses however serious problems in 
general, and in particular if one wants to extend the formalism beyond RPA. First of all any functional is supposed 
to include at least the effect of the short range correlations in the ground state, as produced by the hard core of 
the bare NN interaction. In a microscopic approach this can be achieved by considering a two-body G-matrix in a 
restricted model space as the effective force. For a more phenomenological functional one can include any sort of 
correlation implicitly and have double counting problems. This possible double-counting can be bypassed in a purely 
phenomenological approach, where the force is re-fitted to reproduce the ground state properties and the single particle 
observables after the particle-vibration coupling is introduced at different level of approximation [^ . Notice that in 
the particle-vibration coupling model the particle-particle and hole-hole scattering processes are not considered, since 
the dominant correlation is supposed to be determined by the density fluctuations. 

In the Hartree-Fock approximation the occupation number rii equals 1 for occupied states (holes) and 0 for non 
occupied states (particles). A partial justification of usin g th e HF energies for the single particle orbitals, where the 
nucleons can move, stems from the Koopmans’ theorem [27|, valid for the HF approximation, which states that in 
linear approximation the energy for removing or adding a particle coincides indeed with the HF energy calculated 
in the core. However, adding or removing a particle produces a variation on the Genter of Mass corrections (GMG), 
since they depend on the mass number. The theorem is then violated and the CMC must be properly treated [28j| . 
but in this presentation we are not going to treat this problem. Other physical effects can be relevant for the correct 
description of the single particle potential. If one particle is added or removed, the one-body potential and the single 
particle energies are expected to be modified (’’core polarization”). Our strategy is to incorporate these physical effects 
on the same footing as the dynamical effects on the single particle motion due to the particle-vibration coupling. This 
will be apparent along the development of the formalism. 

We will use the Hamiltonian of Eq. (|3]) to derive the equation of motion of the single particle Green’s function and 
its coupling with the phonons of the nuclear system. To simplify the presentation and the formalism we will consider 
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a generic single particle basis. The normal product can be expanded as follows 

= V'feV'/V'nV'm + < >0 “ < V'feV'm >0 

- < V’/V’n >0 + < V'iVm >0 V'll/’n (6) 

+ < i’ii’rn >0< V'/l/’n >0 “ < V'feV'n >0< V'/V'm >0 

The subscript 0 indicates mean value in the independent particle ground state. It is easy to verify that in the single 
particle basis that diagonalizes the one-body hamiltonian the ground state mean value of the normal product vanishes. 
It will vanish for a generic single particle basis, provided the mean values are still considered in the independent 
particle ground state. As we will see, the formal introduction of the normal product has a relevant consequence in the 
development of the formalism. In fact the single particle self-energy, appearing in the Green’s function beyond the 
HF approximation, introduces modification of the single particle density matrix and of the static mean field acting on 
each particle. These modifications are embodied in a definite term in the equations of motion, see Eqs. (llOllip below. 
In this way one achieves a clear separation between the static polarization of the core, which affects the single particle 
energy and wave function and the dynamics of the nuclear system in the particle-vibration coupling scheme. 

We introduce the single particle Green’s function Gpq{ti,t 2 ) according to the usual definition 

Gpq(ti,t2) = -i < T{il;p(ti)tljl{t2)} > (7) 

where the creation and annihilation operators are in the Heisenberg representation and the mean value is in the 
correlated ground state. The equation of motion for G can be derived following the standard procedure. As shown 
in the Appendix the equation can be derived using the functional derivative method [sol 1^ . In this scheme one 
introduces a single particle external potential 0(t), that depends on the time t, and the dynamical evolution of the 
system under this external potential is analysed. The functional derivatives of different many-body quantities, in 
particular the single particle Green’s function, introduce systematically the correlation functions that characterize 
the dynamics of the system. In particular one has 

- * < > Gnq{ti,t2) ( 8 ) 

0(pkm[tl) 

When all the functional derivatives are calculated at = 0, the equation of motion of the correlation functions of 
the many-body system, in absence of perturbation, are readily obtained. In this way one gets a set of equations that 
couple the correlation functions that is closeT and it can be the starting point of different schemes of approximation. 
The method was previously used in refs. [l3l l32l |33|. In this paper we proceed further in the development of the 
formalism and we will discuss in the next sections its extension to the case of a density dependent force. Details of the 
derivation can be found in the Appendix, while here we describe the main features of the formalism and the resulting 
basic equations. The equation of motion of the single particle Green’s function can be written in the form (Dyson’s 
equation) 

g _ _ 

^ Gpq {tl , ^ 2 ) ^ ^^ ^pp'Gpf q{t\.^t 2 ^ Afp/i (tl, t3)G/ig (ts, ^ 2 ) = <^(^1 i2)Spq (9) 

^ p' h 

where epp> is the the single particle energy matrix, previously introduced, and all the higher order contributions have 
been collected in the self-energy Mph- In this equation, and in the following, a line over a time variable indicates 
integration over that variable (following the notation of ref. 0). The self-energy Mph can be further separated in 
the sum of two terms 


Mphitips) — Uph{ti)6{ti — t^) + Mphitijts) (10) 

The first one is the polarization potential U [H, [s^, 

E < k p\ f \m h >A [< tplitptprnitl) > - < >0 ] (H) 

km 

where the brackets < ... >o indicate mean value in the ground state of Hq, while < ... > indicate mean value in the 
correlated ground state. This polarization potential takes into account the effect of the fluctuations of the density 
matrix with respect to the HF one. Therefore the residual interaction contributes to the static single particle potential 
only as a correction to the original Hartree-Fock single particle potential. This is a direct consequence of the normal 
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product form of the Hamiltonian of Eq. o, which ensures that the residual interaction has vanishing mean value 
in the unperturbed ground state. It is important to remark that the polarization potential is non zero whenever the 
calculation is beyond the original approximation for the single particle energies that appear in the free Hamiltonian 
Hq of Eq. ([S]). Indeed it involves the difference between the higher order density matrix and the one corresponding 
to the one-body part of the Hamiltonian. The correlated density matrix can be evaluated, in principle, for any 
given approximation for the Green’s function. This contribution is just due to the correction to the original static 
single particle potential produced by the correlations in the ground state, as previously anticipated, as well as the 
rearrangement of the single particle potential if one particle is added (or removed) to the nucleus (the static "core 
polarization”). 

The second term in Eq. m describes the dynamical contributions not included in the static single particle 
potential. It can be written 

Mph{ti,t 3 ) = Mp^{ti,t 3 ) + < k p \ f \ h m >a< > S{ti - fa) (12) 

km 


where 




nrij 


(13) 


The self-energy part M' is the one that includes the dynamical polarization processes (dynamical ’’core polarization”). 
In fact the quantity W can be written 


Wpn-tj{ti,t 3 ) =< i p\ f \ j n >A 5(ti - ts) - * < k p \ f \m n >A ^7nk-rs{ti,tb) < r i \ f \ s j >A ( 14 ) 

kmrs 


and H is the particle-hole response function of the nuclear system 


> = - < 4-0 |r{p', (t)^ (t')} I ^0 > 

Prs{t) = < Ipllljs > 


(15) 


It is through the response function H that the phonon degree of freedom is introduced in the formalism. Its coupling 
with the particle degree of freedom (i.e. the Green’ s function) that appears in Eq. (I13|) describes the elementary 
process in which the fluctuations of the single particle potential affects the particle propagation. The vertex function 
r takes into account the effects not included in this coupling process of the particle with the phonon. These effects 
are usually indicated as ’’vertex corrections”. However, the present formalism describes these effects still in terms of 
more complex processes of particle-vibration coupling. As shown in the Appendix, the phonon propagator H and the 
vertex function T satisfy the following two relations that close the set of equations for M, H and T 

nmfc^r's(t,t ) = (t, tl )Ggfc (^2 j tjTpg; (tl , ^2 ! ^3) X 

-t') -iJ2tu < i t \ f \j U >A 


^nh;2j{t3,t4;t5) = SniShjS{t3 — t4)S{t3 — 15 ) — J2rspq Grpib, t8)Gqs{t9, tr)rpq.ij {tg, tg; t^) 

Together with the Dyson Eq. (ED and the ex pres sion m for I^ they form a closed set of coupled equations. Equations 
(I13I14I15I16|1 has been introduced by Hedin IdJ . see ref. for a more recent presentation. In the realm of nuclear 
physics this scheme has been considered in ref. [35|. We stress that the vertex function T is irreducible, i.e. it does 
not contain any polarization propagator starting at the vertex of time In fact this would produce an over-counting 
of contributions in the equation (I16|) for the phonon propagator and unphysical higher order poles 13611. which are in 
fact not present in the expansions of these equations proposed in the literature, see e.g. ref. (^. l37|| . As already 
mentioned, it is in this sense that, in in the terminology of Ref. the vertex function T includes only the 

” vertex corrections”. 

In the integral equation for T the kernel involving the functional derivative of M with respect to G can be considered 
as the effective interaction that determines the vertex corrections. Solving these coupled equations is equivalent to 
the requirement of the self-consistency among M, H and T. In particular one should solve the Dyson equation for 
G. It is important to notice that if the equations are solved in an approximation that goes beyond the independent 
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particle approximation, implicitly some correlations are introduced in the ground state. In general, in the Green’s 
function formalism the explicit ground state wave function is not introduced, but the one-body (from the Green’s 
function) and the two-body (from the phonon propagator) density matrices can be extracted. As mentioned before, 
these correlations modify also the single particle static potential through the polarization potential U of Eq. m- 
The single particle potential will not correspond any more to the one-body part of the Hamiltonian. It is still possible 
to get the energy of the correlated ground state by means of the Koltun’s sum rule [s^, since it involves only the 
single particle Green’s function. However this is not within the scope of the present paper. 

We close this section by mentioning the possible extension of the method to include the pairing correlations. 
Besides the external potential (j) coupled to the density, one can in fact introduce an external source coupled to 
the pairing fields ipitpj and Starting from the anomalous Green’s functions, —i < T{^p(ti)'0q(t2)} > and 

-* < T{V^t(i^)^t(i2)} >, one can generate the propagators and vertex functions pertinent to the pairing correlations, 
and the corresponding extended set of coupled equations. This extended formalism will be discussed elsewhere, and 
we restrict here to the normal systems. 


III. EXPANSION IN DIAGRAMS AND THEIR ESTIMATE 

Once the general formalism is established, one can formulate different approximation schemes. Generally speaking 
one can solve this set of equations by iteration, starting from a guess on the first order approximation. One can start 
from an approximate expression for T and substitute it in Eq. m and (HU, which provides an approximation for 
M and H. The functional derivative of this self-energy is then used to get a new approximation for T. Alternatively, 
one can start with an approximation for M, perform the functional derivative and continue the procedure. One could 
wonder about the second term in the expression (HU for M, which is a mean value not yet specified in the first step of 
the iteration. However, we will see explicitly in the following that this term is cancelled by an identical term arising 
in the expansion of M. This cancellation is exact, and no ambiguity is present in the formalism. 

The iteration provides an expansion of the basic quantities G, H and T, and each term of the expansion can be 
represented by diagrams. At each step of the iteration the so-obtained approximation for T, M and H is conserving, 
i.e. the local current is conserved. As shown in ref. [s^, this property follows from the fact that the effective 
interaction in Eq. (IT6l) for T is obtained by functional derivative of M. Notice that the Eq. (HID for H is also modified 
at each step of the iteration. In principle the starting point is arbitrary, since all possible terms of the expansion are 
automatically generated in any case. 

As suggested by the above discussion, the simplest approximation is to take M = 0 and T = 1. In this way one 
gets the Random Phase Approximation (RPA) for the phonon propagator H, with the single particle Green’s functions 
calculated in the independent particle approximation. This is known to be a conserving approximation. 

In this approximation the set of equations is solved in a consistent way. In fact, in our scheme the single particle 
energies of the one-body Hamiltonian are fixed and their functional derivative vanishes. The RPA describes with 
good accuracy the low lying and high lying collective vibrations. However the corresponding widths require the 
introduction of more complex configurations beyond the standard RPA [s^, see later. One can then continue the 
iteration procedure by inserting P = I in the general expression (USD for the self-energy, perform the functional 
derivative to be inserted in the equation for P, and so on. 

Of course a standard perturbation expansion in terms of the interaction / could be also performed, but the scheme of 
the described iterative procedure involves both the correlated single particle Green’s function and phonon propagator 
(at the given order of the expansion) and at each level of the iteration the corresponding approximation is conserving. 

To develop the PVG scheme within the present formalism a convenient starting ansatz for the vertex function is 
the one that includes explicitly the phonon degrees of freedom at the lowest level. This is depicted in Fig. [I]and its 
analytical expression reads 

= SriShjS{t4 — t5)S{t3 — t^) 

A S/s S/cife2fc3fe4 ^ kir \ f \k2l IIfe2fei;/c3fc4 (G) is) < ^3 s I / I ^4 h > a Gii{t4,t3)Gjs{t5,t3) 

(17) 

To see the connection of this ansatz with the corresponding one for the self-energy, we substitute Eq. (HZD for P in 
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Eq. ([T3l) to get the expression for M, according to Eq. Let us first consider the terms with at most one phonon 

Mph(ti,h) = < r p \ f \ hn >a -t 3 ) + 


+ 1 


2 ^ nr ^-^k\k 2 k^k 4 : 


( < * pI / Ij >A Y.is < kir \ f \ k2l >A n/c2/ci;fc3fe4(t4,i3)x 


X < kz s \ f \ki h >A Gli{ti,ti)Gnr{tl,t4) Gjs{tl,h) + 


(18) 


+ < ki p \ f \ k2 n >A ^k^k^k3k4,{il,t■3) < k:i r \ f \ kih >A 


^ 3 )) + 


+ 5 T.nr < np \ f \hr >A< > S{ti - tg) 

As already anticipated, the first and fourth terms on the right hand side cancel out, and only the two terms with 
one phonon appears in the expression. In the first of these two terms, inside the square brackets a factor appears 
that can be schematically indicated as of the type GGfU. It is then possible to use the general equation (fTHI) for n, 
with r given by the first term only (delta function), which corresponds to the RPA approximation. Schematically the 
relation GGfU = —11 + GG holds true, and explicitly 


^ (^1 ’ (^4? ^ ^ 1^1 f 1^2 ^ ^ A n /,;2 jfcg ^4 (^4, ^ 3 ) — 

— -\- Gnfca (^1 5 ^3)f^fc4Z(^3) ) 


(19) 


Once this expression is inserted in Eq. m, the first term on the right hand side of Eq. m sums up with an identical 
term appearing in that expression, and finally for M one gets 


Mph{tl,h) = J2nrkik2k3ki < kl p\ f \ k2 n >A Uk2kik3ki{tl,t3) < ks r \ f \ k4 h >A Gnr{ti,t3) 

( 20 ) 

~ 2 '^nsijk^ki ^ / 1^7 Gnks (^1; t3)Gkii{t3, ti) < k3 S \ f \ ki h >A Gjs{tl , ^ 3 ) 

The corresponding diagrams are depicted in Fig. The reason why it is necessary to subtract the one-bubble 
diagram from the one-phonon one has been discussed in previous works, see ref. and the review paper Q. If one 
considers that the phonon in the RPA approximation is the sum of many-bubbles diagrams, see Fig. |31 one notices 
that in the one-phonon diagram for the self-energy M, as reported in Fig. [5J the first (one bubble) diagram contains 
two equivalent lines, and therefore the symmetry factor 1/2 is necessary. This is not the case for all the diagrams 
with more than one bubble, where no equivalent lines appear, and therefore the symmetry factor 1 has to be used. 
The subtraction of the one bubble diagram to the one-phonon diagram produces then the correct symmetry factor 
for the first term in the RPA series. 

The explicit expression of the first term in Eq. (1201) . which includes the particle-vibration coupling, can be obtained 
in general by Fourier transforming to the energy representation and expanding both the single particle Green’s function 
and the phonon propagator in their spectral representation. Assuming for simplicity that the Green’s function is 
diagonal, one gets 


= T,n{k} < ki q \ f \ k 2 n >a[ Eac 

<'I'o|Pfc3fc4 |A>< >\<P\lPr,\^o >1 


UJ — U)\ —UJa 


E 


A/3 


bj+ux+ujp—irj 


< k 3 n \ f \k 4 q >A 


( 21 ) 


where A labels the excited vibrational states with energy ui\ of the nucleus with atomic number A, and a, j3 the states 
of the A -I- 1, A — 1 systems with energy uja,uip, respectively. The states a, (3 coincide with the HE particle and hole 
states, and the expression of Eq. (|^ can be written 


Mq(u}) = 


E 

pA 


\FqV 


A 12 


\F^h? 


ui — ui\ — uip + ip 


^ uj + ujx + oJh - ip 


where p and h label particle and hole states, respectively, and 


Fqp = E < ^0 1 A >< fci g I / I p >A 

kik2 


( 22 ) 


( 23 ) 










and analogously for Fhq- 

Since the excited states A are described in the RPA scheme (the starting approximation for our iteration procedure), 
one recovers the standard expression of the diagram, see ref. [llj |. because in this case the matrix elements < 
^olpI.fcJA > can be identified with the amplitudes and [s^ . In ref. [4l| it was shown that taking for a; the 
unperturbed value ujq in the expression for Mq(uj) can be a good approximation. 

As already noticed, two identical terms appear in the expression (jlSIl . each one with a factor 1/2, that sum up to 
the final result ((211)1 . Besides these two terms, one obtains also two-phonon terms for the self-energy. Their graphical 
representations are reported in Fig. (T) Other time-ordering are possible, but we limit the discussion to these two 
cases. The origin and meaning of the diagram (b) can be understood looking at the inner structure of the phonon. In 
fact, if one considers the RPA series that defines the phonon, see Fig. [31 inserted in the diagram of Fig. [51 each bubble 
contains a particle (hole) line that can have the same quantum number as the particle (hole) coupled to the phonon, 
as illustrated in Fig. |4l:, where a bubble of the phonon is put in evidence. To keep the Pauli principle correctly one 
has to consider also the diagram obtained by exchanging these two lines, which is indeed diagramlUD. This diagram is 
therefore demanded by Pauli principle. All the above discussion shows how the expansion is able to take into account 
Pauli principle and to include the correct symmetry factor, which is not trivial once the phonon degrees of freedom 
are explicitly introduced at microscopic level. 

We notice that the choice of Eq. m for the F function modifies the integral equation for the propagator 11. It is 
useful to introduce the polarization propagator P, which is defined in the present formalism as 

Y*rnk;ij = E Grnp (t, ti)Gqk{t2jt)Tpq.ij{ti^t2]ts) (24) 

pq 

It appears as a kernel of the integral Eq. (jl6D for the phonon propagator 11, which in fact can be written 
IImfc,rs (^5 ^ ~ ^ ^ Pmk\ij ^ ^ ^ ^ ^ A ^ut,rs (^ 3 ; O ] 

ij tu 

The polarization propagator corresponding to this choice for the F function is depicted in Fig. O The two diagrams 
(a) and (b) introduce an implementation of the RPA approximation. Notice that the single particle Green’s function 
G has to be calculated with the self-energy just discussed. This is indicated by the thick lines that represent the 
propagator G. In the sequel the self-consistency of the Green’s function will be assumed in any case. Diagram (b) 
includes the coupling of one particle-one hole excitations to two particle-two holes configurations, described as one 
phonon plus a one particle-one hole excitations. In the present scheme the two correlation effects, on the self-energy 
and on the phonon propagator, have to be taken together. This prescription appears as a consequence of the conserving 
structure of the basic equations, as discussed in Section 2. These effects have been considered at perturbative level 
in the study of the Giant Resonances widths in reference [s^, and, using Skyrme forces, in refs. A similar 

scheme has been introduced in ref. [1^, H^, and applied within the Time-Blocking Approximation [47l |. 

It has to be noticed that the structure of the diagrams of Fig. [5] includes the coupling of particle-hole states to two 
particle-two hole states. As such it is physically equivalent to the second RPA [^, [^ for the phonon excitations. In 
the present scheme the equation for the phonon can be written in the RPA form, but with the addition of an energy 
dependent effective force [i^, that comes from the virtual excitation of two particle-two hole states. A detailed 
comparison of this formulation of the second RPA with others [H, is outside the scope of the present work. We 
only notice that the single particle self-energy modifies also the density matrix and consequently, as already stressed, 
the static single particle potential through the one-body interaction U of Eq. ([TTI) . It has also to be noticed that 
these modifications on the phonon can be seen as the introduction of anharmonicity in the vibrational modes due to 
their couplings. 

A second iteration can be performed, which will produce a set of other diagrams, but we will not analyse in detail the 
additional terms generated in this way, since the applications will focus on the first order diagrams. Here we consider 
a general problem that is encountered for the evaluation of the functional derivative of the self-energy. In general the 
derivative ’’cuts” a single particle line within a given diagram for M. This line can appear explicitly in the diagram, 
but it can be also implicitly present in the phonon propagator. It is therefore essential to find a method to perform 
such a derivative ’’inside” the phonon. This can be done by displaying the expansion of the phonon propagator in 
terms of single-particle Green’s functions and interaction matrix elements. Let us illustrate the procedure in the case 
of the RPA approximation. In this case the phonon propagator is the sum of a series of diagrams with increasing 
number of ’’bubbles” , see Fig. [31 The effect of the functional derivative can be figured out by realizing that the Green 
’s function to be ’’cut” can be present at the beginning of each diagram, at the end of it or in the middle. Once the 
cut has been performed, the remaining parts of the diagram can be again summed up, and the phonon propagator 
is again constructed. If the derivative is in the middle of the diagram, on both sides of the cut the phonon can be 
reconstructed. In the other cases this can be done only on one side. Furthermore for the first diagram with only one 
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bubble no phonon can be reconstructed. The functional derivative of the phonon propagator can be then depicted 
graphically as in Fig. [51 The corresponding analytical expression reads 

^ SGraltaM) (^1 ’ ^ 4 )^ = 

= T.pqtu < Pt\f \qu >A + 


+ Y.pqtu ^hk-.pq{tlM) < Pt\ f \qU >A 


<5Grs (^ 65 ^ 7 ) 


(26) 


+ Y.ijab'^]k-ab{tl,h) < ai\ f \hj >A <^^ 1 / 19 ^ >A^ut-lm{t2,U) 

where is the free polarization propagator, and its functional derivative corresponds to the first diagram of Fig. |6l 
This result is in line with the analogous expression derived in reference M- The functional derivative of the phonon 
propagator is essential if the approach has to be kept at a fully microscopic level. 


IV. DENSITY DEPENDENT FORCES 


The diagrammatic expansion introduced in the previous sections implements the formalism of the Nuclear Field 
Theory |5fll - l^ to the case of a realistic effective two-body forces on a microscopic basis. Further developments are 
necessary if modern effective forces or functionals are considered. In general the effective force depends on the density 
or the density matrix and the formalism must be generalized to this case. In minimizing the EDF one has to consider 
the functional variation of the force, and the mean field (Hartree-Fock) potential V does not involve only the effective 
force but also its functional derivative 


< fc I V I fc' > 


^^Pkk’ >0 


1 

2 


'Zllj 


pI + pU 


j J lk;jk' 




(27) 


Solving self-consistently the resulting Hartree-Fock equations with this effective single particle potential is equivalent 
to minimize the full functional of Eq. O- Notice that the first two terms in square brackets are equal. 

The dynamics of the density fluctuations is determined by the energy surface around the minimum. If the fluctua¬ 
tions are not too large, one can expand the functional to second order in the fluctuations 5pij 


E — Eq + \ 7 -) 5pij + — ( 7 - 7 -) 

^ ^ Pij ^ ^ ^ ^ PsfV Pnv ' 


S^V 


^Pst^Pu 


Spst^Pu 


(28) 


Also in calculating the second derivative one has to take into account the density matrix dependence of the effective 
force 


< U S \ K \ v t >A = ( 


s^v \ 

SpstSpuv ' 0 


= M ft,., + ft,. + [E« ( 




fA. 

^0 ( ^ hv,k 


(29) 


The quantity k can be interpreted as the interaction between two density fluctuations. This is a well known result, 
usually considered in connection with the standard RPA approximation for the vibrational states [38l| , and we report it 
here for later discussion. As such Eq. (1291) gives the matrix elements of the basic interaction between two elementary 
particle-hole excitations. If the effective force is density independent, only the first two terms in the curly brackets 
are present and the matrix element, as expected, is correctly anti-symmetrized and it coincides with the previous 
two-body matrix element in the case of density independent forces (notice that the first two terms are actually equal). 
The density dependence introduces the additional ’’rearrangements” terms involving functional derivatives. These 
terms are not automatically anti-symmetrized, despite we started with the anti-symmetrized matrix elements of Eq. 
©. This is due to the already mentioned lack of anti-symmetry with respect to orbitals that appear in the density 
dependence of the effective force. In principle one should then anti-symmetrize these matrix elements. This is a 
general problem for the density dependent effective forces. Relying on the first line of Eqs. (071) and (0^ . we shall 
develop the formalism for a generic energy density functional. In that case, we are not supposed to trace back any 
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direct or exchange term of any starting interaction V. Of course, if such interaction exists like in the case of the Gogny 
force and most of the Skyrme sets, the equations on the second line of Eqs. (HZl) or (1^ can be of some help. The 
identification of the residual interaction with the second functional derivative can be further clarified if the density 
dependence of the functional is considered explicitly. The single particle density p(r) can be written in terms of the 
density matrix (neglecting spin-isospin labels) 


p(r) =< V'(r)'^ lAW >= X! < 

ij ij 


(30) 


For a fixed single particle basis {(^}, the variation of the density matrix implies a definite variation of the density. It 
follows that the functional derivative of V with respect to the density matrix can be written 


^ Pi j 

Similarly one gets 


= / d-^r 


9p{r) 

dp{r) dpij 


= / 


(5F 

( 5 p(r) 


=< i 


SV 

6p{r) 


j > = <i \ V \ j > 


S^V 


dpij dpim 


^ I I Sp{r)Ip{r') ^"^^'** <?i'm(r') =< i l\ k \ j m > 


(31) 


(32) 


This is the general form of the effective density-density interaction associated with a given EDF if one can stick to the 
harmonic approximation, i.e., stop the expansion of the potential with respect to the density at second order. This 
can be justified in specific cases by looking at the transition density associated with the vibrational states, which can 
be small with respect to the corresponding unperturbed density, see e.g. figure 1 of ref. (s^j . 

The harmonic approximation, i.e., the truncation of the expansion at the second order, may be inappropriate in 
some cases; however, one could in the same way as in Eqs. (j31l32|l identify the matrix elements associated with third, 
fourth, etc. order variations of the potential. 

At this point the development of the formalism can follow exactly the same lines as in the case of a density 
independent force. We adopt the same interaction that describes the density-density coupling also for the particle- 
vibration coupling. This goes along the line of Ref. where it is shown that the matrix elements for the 3 particle-1 
hole coupling processes have indeed the form of Eq. (l32ll (in all cases, except for the quantum numbers J'" = 0+ of 
the excited particle-hole). The problem of double counting, discussed previously, arises also in this case, and again a 
refitting of the force parameters should be envisaged. 

However, an additional caveat has to be considered. Since the equations of motion are derived using the form of Eq. 
(|29|) for the matrix elements of the two-body interaction, the processes that determine the single particle dynamics 
and the nuclear vibrations can include particle-particle scattering, for which the density-density interaction matrix 
elements of Eq. (1291) does not look appropriate. In these processes, in fact, the particle-particle or hole-hole interaction 
does not involve density fluctuations . The hamiltonian of Eq. o necessarily produces these elementary processes. 
This feature is intrinsic to the hamiltonian formalism, where the two-body matrix elements are fixed from the start 
and they cannot depend on the processes where they are involved. 

To see how particle-particle scattering arises in the present formalism, let us consider the functional derivative of 
the second diagram (the ’’bubble”) of Fig. [21 which is performed at the second iteration of the expansion procedure. 
If the functional derivative acts on the hole line of the diagram, the vertex function develops the diagram (a) of Fig. 
0, which in turn generates the diagram (b) for the self-energy. The interaction at the center of this diagram is indeed 
a particle-particle scattering, which is expected to be described by a matrix element different from (I32|) . 

The particle-vibration model assumes that this processes have negligible effects, since the particle-vibration coupling 
and the particle-hole matrix elements are considered the dominant ones. In this case no conceptual difficulty can arise. 
Along the same talking, the approximations obtained in the expansion are conserving to the extent that these processes 
can be neglected. 

In any case we will restrict the applications to the lowest order approximation, where this problem is absent, and 
we will leave a throughout analysis of this point to future works. 


V. APPLICATION TO ‘^°Ca. 

The lowest order dynamical contribution to the single particle self-energy is given by the two diagrams of Fig. [2] 
Notice that they have to be subtracted. As discussed in the derivation of this result in Sect. Ill A, the subtraction of 
the bubble diagram is needed once the phonon degrees of freedom is introduced in order to fix the correct symmetry 
factor. We will study the relevance of this additional diagram in the specific case of the nucleus ^°Ca. It should be 
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TABLE I: Single particle energy levels of for the Skyrme force SV. The column HF contains the Hartree-Fock energies, 
Eph and Ehub indicate the contribution of the one-phonon diagram and the bubble diagram of Fig. [2l respectively, Ecor their 
difference. The column labeled E contains the corrected level energies and E^xp the experimental ones. Energies are in MeV. 



HF 

^ph 

Ebub 

Ecor 

E 

Eexp 

ld5/2 

-27.67 

-0.76 

0.53 

-1.28 

-28.95 

-22.39 

2sl/2 

-19.34 

-4.25 

-0.42 

-3.83 

-23.17 

-18.19 

ld3/2 

-18.95 

-1.58 

0.32 

-1.90 

-20.85 

-15.64 

lf7/2 

-7.54 

-0.87 

-0.30 

-0.57 

-8.11 

-8.62 

2p3/2 

-2.02 

-2.28 

-0.75 

-1.53 

-3.55 

-6.76 

2pl/2 

-0.37 

-1.74 

-0.60 

-1.14 

-1.51 

-4.76 

lf5/2 

1.62 

-1.21 

-0.42 

-0.79 

0.83 

-3.38 


TABLE II: Single particle energy levels of 

for the Skyrme force Sly5. 

Notation as in Table 1. 



HF 

^ph 

Ebub 

Ecor 

E 

Eexp 

ld5/2 

-22.09 

-1.30 

-0.88 

-0.42 

-22.51 

-22.39 

2sl/2 

-17.27 

-3.54 

-2.29 

-1.25 

-18.52 

-18.19 

ld3/2 

-15.19 

-1.65 

-0.98 

-0.67 

-15.86 

-15.64 

lf7/2 

-9.67 

-2.34 

-1.50 

-0.84 

-10.51 

-8.62 

2p3/2 

-5.29 

-4.70 

-2.71 

-1.99 

-7.28 

-6.76 

2pl/2 

-3.10 

-4.64 

-2.72 

-1.92 

-5.02 

-4.76 

lf5/2 

-1.28 

-2.71 

-1.12 

-1.59 

-2.87 

-3.38 


noted that this additional diagram had been considered in the past, but not in any of the recent self-consistent PVC 
calculations. Our study, therefore, is the first of such kind with the consistent use of the whole Skyrme central force. 

Let us first consider the case of the density independent Skyrme force SV [l^. In Table I we report the single 
particle energy levels around the Fermi energy. The Hartree-Fock (HF) energies are displayed, together with the 
phonon diagram contribution Eph and the bubble diagram contribution Ehub- The total correction Ecor = Eph — Ehub 
is added to the HF energy to obtain the “dressed” single particle energies E. The calculations have been performed by 
simply setting an upper cut-off of the single-particle energies at 60 MeV: the RPA calculations and PVC calculations 
are done within this model space and no additional selection has been made, so that in particular all RPA phonons with 
multipolarity O'*", 1“, 2'^, 3“ and 4+ are retained. In this way, the phonon contributions are corrected consistently by 
the corresponding bubble diagrams. The self-energy corrections due to the PVC coupling turn out to be all negative, 
both for particle and hole states. The reason of this particular feature is discussed in ref. m- One can notice that 
the correction due to the bubble diagram is not negligible and in some cases it is substantial. 

To study the dependence of the results on the adopted effective force, we consider the more modern Skyrme force 
Sly5 [13 . Since the force is density dependent, we follow the scheme described in Sec. IV. The HF energies and the 
self-energy contributions are reported in Table H. Even in this case, one can notice that the bubble diagram is in 
general a substantial fraction of the phonon diagram contribution and its overall size is even larger than for the SV 
Skyrme force. Therefore, we conclude that the inclusion of the bubble diagram certainly requires careful consideration 
in all future calculations performed in the PVC scheme. 

We do not wish to discuss in detail the comparison with experimental data, let alone to claim that the inclusion 
of the PVC coupling is improving such comparison. We could note, on the one hand, that the r.m.s. devation with 
respect to experiment of the HF energies is 1.28 MeV, and that this value increases to 2.40 MeV with the inclusion 
of the one phonon diagram contribution Eph but it decreases to 0.78 MeV with the addition of the bubble diagram 
contribution Ehub- On the other hand, one must be aware that the results obtained with a zero-range force are very 
sensitive to the choice of the model space and can diverge if no cut-off is set. While we have started to tackle this 
problem [H, H^l, the present results are only meant to illustrate that the lowest-order approximation to PVC for 
single-particle states must include the bubble diagram. 


VI. CONCLUSIONS. 

We have presented a many-body theory of nuclear structure within the particle-vibration coupling scheme, where 
the degrees of freedom are single particles and phonons, each one described by the corresponding propagator. The 
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formalism is developed for a generic effective interaction or Energy Density Functional, which are supposed to be 
defined within a model space. Both density independent and density dependent effective interactions have been 
considered. The theory is devised to go beyond the mean field approximation to all orders. Since the effective 
interaction is expected to include part of the correlations, in principle a refitting of the effective force would be 
mandatory. We have not approached this problem, but instead we have focused on the correct treatment of the 
particle-vibration coupling when the phonon degrees of freedom are treated at a completely microscopic level. The 
method is based on the equations that couple the single particle self-energy, the vertex function and the phonon 
propagator, and are derived within the Hamiltonian formalism. The solution of the equations is by iteration. At 
each iteration one gets an approximation that is conserving Sim , i.e. the vertex function conserves the scalar 
current and the total momentum. The terms so obtained by iteration can be represented by diagrams. In this way 
both static and dynamical corrections beyond the mean field are systematically introduced. We must stress that the 
diagrammatic expansion takes into account systematically both the Pauli principle and the possible over-counting of 
contributions. We have shown that with a simple one-phonon ansatz for the vertex function one obtains in a natural 
way the lowest order approximation for the self-energy in terms of Fig. [2l 

The general expansion provides a formal extension of the Nuclear Field Theory to realistic effective inter¬ 

actions and it is a non-perturbative solution of the many-body problem based on an effective interaction in a model 
space. Truncating the expansion at a given order provides a conserving approximation within the particle-vibration 
coupling scheme. 

We have presented applications of the theory to the nucleus '^°Ca. Both the density independent Skyrme force SV 
and the density dependent Sly5 were considered and compared. Already at this level the results indicate that it is 
essential to take into account correctly the microscopic structure of the phonon, which implies the introduction of the 
bubble diagram of Fig. [5] to fix the correct symmetry factor for the underlying RPA diagrams that define the phonon 
propagator. The bubble diagram turns out to be a substantial correction to the one-phonon diagram, and demands 
a careful treatment in any application of the PVC coupling scheme. 

In summary, we claim that even at the lowest order a sound PVC calculation should include at least the two 
diagrams for the self-energy in Fig. 2. Although this result has been known for some time, we are not aware of its 
implementation in recent self-consistent calculations of the self-energy with microscopic energy functionals (although 
such correction has been introduced in the total energy calculations of Ref. [13). Our numerical results indicate the 
importance of the bubble correction. More importantly, we pave the way to calculations that go beyond this lowest 
order approximation. 


Appendix A: Derivation of the basic equations 

We give here some detail on the derivation of the set of the three equations that are the basis of the formalism 
described in Section In] In general the equation of motion for G can be written 

— ^(ti t2)Spq y~(„/ epp'G*D'g(tli ^2) 

- *5 Efcmn [< kp\ K\mn >A< > (Al) 

-I- < kp \ K\nm >A< >0 Gkqitl,t2)] +J2l4>plitl)Glqitl,t2) 

where the label A indicates anti-symmetrized matrix element. The two-body interaction term (second line) is directly 
related to the functional derivative of the Green’s function with respect to the external potential (j){t) 

- i < T{l/'l(tl)'i/’m(G)V'»(D)V’J(^2)} > = > Gnq{tut2) (A2) 

It has to be noticed that the mean-value at the right hand side is time dependent due to the external interaction (j)(t) 
and must be kept distinct from the static mean-values that appear in the expansion of the normal product of Eq. (jb]) . 
The former is just the time dependent one-body density matrix of the system. Let us introduce the inverse Green’s 
function G~^ 


'y {ti,t2)Grq{t2,tl) — Gpr (tl, (^2 , G) — SpqS{ti t^) 


(A3) 
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where a line over a time variable indicates integration over that variable following the notation already introduced in 
the text. From Eq. (TMI) one gets 


SGnq{ti,t2) /'+ ^ ri-^\ 


^4^km (^1) 


(tl ) 


(A4) 


If we insert this expression in Eqs. (IA2IA1I) . the equation of motion (lAlIl can be written according to Eqs. (I MT^ . 
where 


Mph{ti,t4) = -4 < k p\ K \m n > Gnr(ti,ta) fe^4) 

^ 9 ^ d(pkm{tl) 


(A5) 


kmnr 


By comparing Eqs. (lAip and (IA4I1 one gets the expression of the inverse Green ’s function G ^ in terms of the 
self-energy 


“ ^ 2 ) - epqS{ti - t2) 

- [M^g{ti,t2) + ^ J2km < k p \ K.\qm >A< > S{ti - t2)] 

- Upq{tl)6{tl - t2) 


(A6) 


where the potential U is given by Eq. m and the quantity in square parenthesis is just M of Eq. m- 

The phonon degree of freedom can be introduced by noticing that the functional derivative of the potential U 
is related to the density-density response function, that can be identified with the phonon propagator. In fact the 
functional derivative of U can be computed according to the general method as in Eq. (1X21) 


= SprShsS{t - t') - * V' < k p \ K\mh ^8r-km{t, t') 

where 11 is defined according to Eqs. (ITSl) . In fact, according to the general result of Eq. (IA2I) . one has 

= * < T{tpl{ti)tl;mitl)lpl{t2)lpnit2)} > - i< V'l (^1 (G ) >< V’J )'*/'« ) > 


(A7) 


(AS) 


— ^14m/c;gn(^l; ^ 2 ) 


The particle-vibration coupling can be introduced using the chain property of the functional derivative, since the 
potential 17 is a functional of the external perturbation cj). Indeed the M' self-energy can be expressed in terms of the 
phonon degree of freedom 


M. 


phiGM) = -i\T.kmnr^j< k p\ K \m U > A Gr,r{h, t4) 

— 2 ^^^nrij ^^pn;ij {^ 1 : i 5 ')Gnr{il: rh',ij {^ 3 : ^ 4 t ^ 5 ') 


(A9) 


where 


Wpn■,^Jitl,h) = Sfcm < k p\ K\mn > {) 




(AlO) 


Together with Eq. dm, this gives Eqs. (I13I14I) . The equations (fTHll for 11 and T can be derived using again the chain 
property of the functional derivative. For the phonon propagator one gets 

TT (M M \ _ SGnq{t2,t2 ^ ) 

^^mk-qn(l'l, l'2) — S4,km(ti) 


_ SGriq(t2,t2 y SUjj (< 3 ) 

“ ^^3 SUij(i^) ^ HkmiG) 

~ '^ijrs ^^r{ti,t4)Geq{t4,t^)Trs;ij{t4,t5',t2~^) X 


(All) 
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where the definition of F in Eqs. (lAlOl) and Eq. (IA4I) have been used. When Eq. (IA7|) is substituted in Eq. (lAllI) 
one gets the first one of the constitutive Eqs. (fTBl) . In Eq. (lAllI) t 2 = h + Vi where 77 is a positive infinitesimal, 
to ensure the correct order of the creation and annihilation operators in the equal times Green’s function. For the 
vertex function F one gets from Eq. (IMl) 

F„?i;jj(tl, ^2^3) = “( Slhlits) 






h) - ( 




(A12) 


= 5r^5hJ5{tl - h)5{t2 - h) + llrspq{^i§^^^)Grp{tA,tQ)Gqs{t7,h){^^ 


(t6.t7)^ 




and with the definition for F in Eqs. (lAlOll one obtains the second of the constitutive Eqs. (ITH)) . This completes the 
derivation of the three basic equations of the formalism. 
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FIG. 1: Diagrams representing the starting guess for the vertex function in the iteration procedure described in the text. 
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FIG. 2: First order diagrams for the self-energy beyond HF, consistently with the ansatz for F in Fig. 1. Notice the subtraction 
of the ’’bubble diagram” : its symmetry factor and other details are discussed in the text. 
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FIG. 3: The diagram series for the RPA approximation. 
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FIG. 5: Polarization diagrams as obtained at the second iteration step. The thick lines indicate that the Green’s function are 
considered consistently beyond the HF approximation. 



21 



FIG. 6: Pictorial description of the effect of the functional derivative on the phonon propagator. The first picture on the 
left corresponds to the functional derivative of 11°, which appears at the left hand side of Eq. ([26]). See the text for further 
explanation. 
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(a) 


(b) 




FIG. 7: (a) Diagrams for the vertex function generated by the functional derivative of the second diagram of Fig. @.(b)The 
corresponding self-energy diagram, where the interaction in the middle corresponds to a particle-particle scattering process, 
see the text. 







